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1 Introduction 



There has been a long history to study an asymptotic solution of Hamilto- 
nian flows by means of singular perturbation methods such as the averaging 
method and the method of multiple time-scales. A Hamiltonian flow can be 
reduced to a symplectic discrete map called the Poincare map, which has the 
lower dimension than the original flow and is, therefore, extensively studied 
PQ[2]. However, neither the averaging method nor the method of multiple 
time-scales may be immediately applicable to symplectic maps. 

The perturbative renormalization group (RG) method developed recently 
may be a useful tool to tackle asymptotic behaviours of discrete maps as well 
as flows. The original RG method is an asymptotic singular perturbation 
technique developed for diffential equations [3]. Secular or divergent terms 
of perturbation solutions of differential equations are removed by renormal- 
izing integral constants of the lowest order solution. The RG method is 
reformulated on the basis of a naive renormalization transformation and the 
Lie group [3]. This reformulated RG method based on the Lie group is easy 
to apply to discrete systems, by which asymptotic expansions of unstable 
manifolds of some chaotic discrete systems are obtained [5]. The extension 
of the RG method to discrete symplectic systems is not trivial because the 
symplectic struture is easily broken in naive RG equations (maps) as shown 
in Ref . (H] , while the application of the RG method to Hamiltonian flows does 
not cause such a problem as the broken symplectic symmetry except a special 
case [7j. 

The application of the RG method to some non-symplectic discrete sys- 
tems has been attempted in the framework of the envelope method How- 
evever, the method, if applied to a symplectic map, would give only a naive 
RG map which breaks the symplectic symmetry. To recover the symmetry 
breaking, not only the paper [S] has introduced a process of the "exponentia- 
tion" which gives us analytical expressions of some physical quantities jHUH], 
but also another paper has used a symplectic integrator by taking the 
time-continuous limit of the time-step parameter. Furthermore, S.I. Tzenov 
et al have showed that the exponentiation procedure can be successfully ap- 
plied in the Henon map|llj. However this procedure may seem somewhat 
artifical, and resonant islands have never been studied by the RG method. 

The main purpose of the present paper is to present a general RG pro- 
cedure to preserve the symplectic structure in RG maps and to obtain cor- 
rectly reduced symplectic RG maps. In this paper, this procedure is called a 
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symplecticity-preserving RG procedure, which consists of the following two 
steps. First, using the reformulated RG method jlj, we get a naive RG map 
near elliptic of some symplectic discrete systems. The naive RG map pre- 
serves the symplectic symmetry only approximately and fails to describe a 
long-time behavior of the original maps. Second, in order to recover the 
symplectic symmetry by means of taking account of higher orders in small 
parameters in the naive perturbation result. Then we equate the naive RG 
map to the appropriately discretized Hamiltonian flow, which is done by 
using the Liouville operator. On this procedure we identify the small param- 
eter e with the time-step, which yields a canonical equation in the limit of 
e — > 0. This process and a symplectic RG map are called, respectively, a 
symplecticity-preserving procedure and a symplecticity-preserving RG map. 
A reduction of the Liouville equation in a non-integrable Hamiltonian flow 
has been developed in Ref . [T2"] . 

In section 2, a long-time behaviour of a simple linear map is analyzed 
in order to elucidate both the broken symplectic symmetry in a naive RG 
map and our regularization process. In section 3 , a symplecticity-preserving 
RG map is obtained near elliptic fixed points of a two-dimensional nonlinear 
symplectic map even when a resonant island chain appears. In section 4, we 
mention the advantages of our newly developed RG method as the summary 
of this paper. 



2 Linear symplectic map 

It may be instructive to analyze a linear symplectic map, which is exactly 
solvable, (x n ,y n ) i— > (x n+1 ,y n+1 ): 

x n+1 = x n + y n+1 

y n+1 = y n -ax n + 2e.Jx n . 

This map can be rewritten as 

L e x n = e2Jx n , (1) 

where e is the small parameter, J is the real parameter, 9 and Lgx n are 
defined as 

costf = (l-a/2), (2) 

Lgx 11 = x n+1 -2x n cos6 + x n -\ (3) 
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respectively. We assume that the map (JTJ) has an elliptic fixed point at the 
origin (0,0). The linear map (JIJ) has the following exact solution x 1 ^ : 



x E 



A exp i arccos(cos 9 + eJ)n I + c.c 



Aexp 



i\e + e 



-J 



sin 9 



+ 6* 



cos 9 ( J \ 2 



2sin# \sin# 



n 



+ c.c, (4) 



where v4(g C) is the complex "integration" constant and c.c. stands for the 
complex conjugate to the preceding terms. 

Let us derive an asymptotic solution of the map for small e by means 
of the RG method. Substituting the expansion 



x n = x (0)n +£a; (l)n +£ 2 ;r (2)n + 



(5) 



into Eq.(£Q), we have 

L e x^ n = 0, L e x^ n = 2Jx^ n , L 6 x^ n = 2Jx«" • • • . 



and 



x 



(0)n 



X 



(l)n 



X 



(2)n 



Aexp(i9n) + c.c, 

— nexp(t9n) + c.c, 

sin 9 

-J 2 A 



2 sin^ 9 



2 .cos 9 

n +t -n ) expnOn) + c.c. 

sin 9 ' 



where A(e C) is the integration constant. To remove secular terms (oc n, n 2 ), 
we introduce the renormalization transformation A i— ► A n 

-iJA -J 2 A 



A n = A + e- 



sin 9 



-n + e z 



2 sin 2 



?2 + % 



.cost 
sin 9 



-n 



0{e" 



(6) 



A discrete version of the RG equation is just the first order difference equation 
of A n , whose local solution is given by Eq.flBJ). From Eq.©, we have 



A n+1 — A n 



J 



J 2 



is 



— e 



sin 9 2 sin 2 9 



In + 1 + i- 



.cos 9 



sin i 



\A + 0(e 3 ), (7) 



where A should be expressed in terms of A n . This is done by taking the 
inversion of the renormalization transformation © iteratively 



.4 



l + t e^- + 0(e 2 ))A^ 
sin 9 



(8) 



Substituting (JSJ) into (J7J), we obtain the following RG equation (RG map) up 
to 0{e 2 ) 



A n+1 

of which solution is 

A n -- 



—ieJ 1 

1 + — ~W + 



sin 9 2! V sin 9 



-ieJ 



IE 



, J 2 COS " 
2 sin 3 



A n + C(e 3 ), 



(9) 



— zeJ 1 
1 + — ~w + ^ 



sin9 2! V sin# 



-ieJ 



IE 



, J 2 COS 

2 sin 3 9 



+ 0(e 3 )) A . (10) 



On the other hand, from Eq. (HI), we have A n exactly as 



A n = A°exp 



Z £ 



-J 

sin # 



cos^ / J \ 2 



2 sin 6 1 V sin 9 



+ 



n 



(11) 



Notice that \A n \ 2 is an exact constant of motion while it is merely an ap- 
proximate conserved quantity of the (trancated) RG map (jUJ)- The symplec- 
tic structure is also not exactly preserved in the RG map, that is, for the 
truncated RG map Q up to 0(E k ), we have 



dA n+1 A dA 



* n+l 



dA n AdA* n = 0{e 



(M-i^ 



where k = 1,2, • • •. A* is complex conjugate to A n and should also be a 
canonical conjugate to A n [14 j. Although this fault of the RG map vanishes 
in the limit of e — > 0. 

In order to remedy a fault like this, we take advantage of a crucial obser- 
vation that Hamiltonian flows satisfy the following relation: 



// 



Z(t + {i) = [l + fiC H + + )Z(t) = exp(^£ H )Z(t) 



(12) 



where t is the time variable, [l is a real number, if is a Hamiltonian, Z is a 
canonical variable ((ft, g 2 , •■■) Qn,Pi,Pz, ■■■,Pn) an d 



A' 



7^ \dqj dpj dpj ihij 



(13) 



{{Z,H},H}, 
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where Ch is the Liouville operator defined by the Poisson bracket. We can 
identify the relation (|12j) as the map by defining Z n+1 = Z(t + /i), and 
Z n = Z(t): 

Z n + i = ^ (Z n ;/x)> m^n) = exp(^£ H )Z(t)\ zmzn . (14) 

This is the symplectic mapping associted with the Hamiltonian system. Iden- 
tifying fi — e in (jl4}) . we can obtain both the Hamiltonian system and the 
associated symplectic mapping in the following steps: 

1 Exapanding a Hamiltonian H in (jl2J) in powers of e, H = +eH^ + 
■ ■ •, we obtain the following relation 

Z{t + e) = (l + eC H + £ -C 2 H + 0(s^Z(t) 

= {l + eC HW + e 2 {^ + £ H(2) ) + 0{e i )]z{t). (15) 

2 can firstly be found by taking the limit (A n+1 - A n )/e -> dA/dt, 
(e — > 0) in the naive RG map by comparing the naive RG map with 
()15|). Using and eqauting the naive RG map with the Liou- 
ville operator relation (J15)) . we obtain the appropriate Hamiltonian 
H = + eH^. Similarly, we can obtain the higher Hamiltonian 
H( 3 \ order by order in e. This procedure yields the approxi- 
mate Hamiltonian. 

3 To obtain the symplecticity-preserving RG mapping, we discritize the 
continuous-time system obtained in the 2nd step. The time-step should 
be chosen to be e. 

In the present case, we take the limit e — > in Q, which gives 

dA iJ OH dA* OH 

~~r~ = — : — 7.A = - - = jCij(i)A, — : — = — - = C-u{i)A . 

dt sin# dA* H ' dt dA* H 

Here ^ is from (A n+l — A n )/e, and so on. Then we have 

ffW = -i^, and Cl w A = {{A,H^},H^} = ^. 

sin V sin 6* 
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At this stage, the Liouville operator relation is 



1 + eC RW + e 



r 2 



C 



.4 



A + e 



-UA 



-J 2 A dH&\ 
+ 



sin9 ' " V2sin 2 # ^ &4* Y 
Equating this to the right hand side of the naive RG map (JHJ), we have 



-iJ 2 cos6»|v4| 2 
2 sin 3 9 



The Hamiltonian H = H^ 1 ' + leads to the following equation: 

dH 



dA -U -i J 2 cos 6 A _dH 
~dl ~ s^m9 + 6 2sin 3 # ~ ~dA ' 
The solution of the continuous-time system is 



dA 
~dt 



dA* 



i / —J —J 2 cos 9 . 
A(t) = A O exp z — + g 1/ 
\sint/ 2 sin 6^ 



which yields the symplecticity-preserving RG map: 



.4 



ra+l 



A n exp <^ is 



-J 



-J 2 cos 9 



^sin9 2sin 3 # 

Here the time-step is chosen to be e. Then we recover the exact solution of 
CEJ, is pijl . in terms of the original variable x n . 

Finally, the present method can be summarized in the following diagram: 

A symplecticity-preserving RG method for symplectic maps 



[discrete-time 



[cont inuous-t ime] 



Symp. maps 

naive RG I 

naive RG maps 
(dissipative maps) 

Symp. -Pres. RG maps 



Liouville operator 
> 

Symp.-Pres. discretizations 



RG Eqs. 
canonical Eqs. 
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3 Two-dimensional Nonlinear Symplectic Map 

3.1 Non- resonant case 

Let us analyze the weakly nonlinear symplectic map (x n , y n ) t— > (x n+1 , y n+1 ): 

x n+l = x n + y n+l 

y n+1 = y n - ax n + 2eJ{x n f ) 

or 

L e x n = e2J(x n ) 3 , (16) 

where e is the small parameter, J is the real parameter, 5 and Lgx n are 
defined in (J2|) and (J3J) respectively. Expanding x n as a power series of e 

x n = x (0)n + £x (l)n + x (2)n + 

we have 

L e x {0)n = 0, L e x {l)n = 2J(x (0) ") 3 , L e x (2) " = 6J(x (0) ") V 1)ri . (17) 
and the solution to the perturbed equations to 0(e 2 ) are given by 

x (0)n = Ae i6n + cc ( lg ) 

* (1)n = ~ 3i|A ' MJ ne^ + fl e^ + c.c. (19) 

sm & cos 35 — cos 5 

(2)n f-9 J 2 |A| 4 A 2 , J 2 |A| 4 A / 3 9cos5 \ 1 i0n 

X ~ IT sin 2 5 ™ sin5 Uos35-cos5 + 2^n^J n J e 

(_9iJWA^ J 2 |A| 2 ^ f 12 _ sin 3(9 U ^ 

I (cos 35- cos 0) sin^ 2(cos35 - cos5) 2 I sin5 J J 

+ I 3J ^ 5 \ e 5ie " + c c (20) 

I (cos 55 - cos 5) (cos 35 - cos 5) J 1 J 

Here A(G C) is the integration constant. To avoid resonant secular terms, 
we assume in addition that 

cos 9 cos 35, cos 5 ^ cos 55. 

The construction of the reduced map in a case of near resonance is considered 
in subsection 3.2. In order to remove the secular terms in the coefficient 
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of the fundermental harmonic (exp(zftn)), we introduce the renormalization 
transformation A i— > A n : 

A n = A + e 1 t-^ n 



sin ft 
-9 J 2 \A\ 4 A 

~2 sin 2 ft 
.J 2 \A\ A A{ 



-n 



—i- 



sinft V cos 39 — cos 9 2 sur ft 



+ 



9 cos ft 



n 



(21) 



Following the same procedure as that in the preceding section, we derive the 
naive RG map from (J2*T|) : 



.4 



n+l 



A" + e -^\A n \ 2 A n 
sin ft 



-3zJ, 



9i cos ft 



J 



2! V sin ft 
3iJ 2 



\A 



n\2 



A n 



2 sin 3 ft sin ft (cos 3ft — cost 



\A n \ A A n 



(22) 



which breaks the symplectic symmetry. This RG map should be recovered 
the symplectic symmetry. The naive RG map can be written in terms of the 
real variables A n = A™ + iA%, 



A\ 



n+l 



Q 1 

A1 + e— —{Af + Af)A™ 
sin ft 



+e 
+ 



-If 3J 

^Tlsin 
9 cos ft 

+ 



3 7 1 2 

in ft J 



2 sin ft sin ft (cos 3ft — cos i 
—3 7 

A" + e— — (Af + Af )A n l 
sin ft 



^J 2 {Af + AffA^ 



(23) 



1 

r 

9 cos ft 



+ 



i4» 



The (e 



. 2 sin 3 ft ' sin ft (cos 3ft — cos i 
0)-limit yields 



J 2 (Af + AffA^ 



(24) 



~dT 

dA 2 
dt 



3J 
sin ft 
3J 



sin i 



Ai + A£)A 2 = 
Al + AjjA, 

9 



<9# 



£tf(l)-<4lj 



OH 



where 

H{1) = 3J(Aj + Alf 
4sin# 

This Hamiltonian generates the relations: 

3J 



A l+ e^-(Al + A 2 2 )A 2 



12 

3J 



A 2 + e—-(Al + A 2 2 )A l 
smt) 



According to the general procedure, can be obtained by equating them 
to Ij23jl -(|2 ^ . That is, 



R(2) = f 9cos0 + 3 | J^4| + A§g 

1 2 sin # 3 sin#(cos3# — cos#) J 6 

Then we obtain the approximate Hamiltonian if = + eif ^ 

# = a{Al + Al) 2 + /3{Al + A 2 2 ) 3 . 

3J r9cosfl 3 1 J 2 

" = 4sinfl' ^ = £ l2sin 3 ^ + sin#(cos3fl - cosfl) J IT 

The canonical transformation dA\ A dA 2 = d® A dl: 

A x = y/2I sin 6, A 2 = y/2I cos 6, 
gives us the simple canonical equations: 

ore a/ 

dl _ _ 

¥ ~ ~dQ~°' 
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In addition to this expression, there exists the exponential form 

A = A 1 + iA 2 

2/(0)(sine + icos9) 



2/(0)i(cos9 -zsin9) 



2/(0) exp - i 8a/(0) + 24/3/(0)^ )t - £0(0) + ztt/2 



Therefore, the solution of (e — > 0)-system which is written in the form of 
exponential function is 



A(t) = A(0)exp\ - it 
+24e 



\A(0)\* 
4sin# 2 
9cos# 3 



2 sin sin 9 (cos 36 1 — cos i 



j 2 \a(o) r- 



A(0) exp | - it 



3J 
sin # 



L4(0)p 



9cos# 



2 sin sin 6 1 (cos 39 — cost 



where we use the relation J 21 (t) = \A(t)\ = const. To construct the reduced 
map, we discretize the Hamiltonian flow. Identifying 

A n+1 = A{t + e), A n = A(t), 

yields a symplecticity-preserving RG map 

-3J\A n \ 2 



A 



n+l 



A n exp 



is 



+eJ 2 \A n \ A 



sin# 
9cos# 



2 sin 9 sin#(cos30 — cosf 



(27) 



This map can also be obtained by the exponentiation RG method [§]. 



3.2 Resonant case 

Let us consider a case of resonance in the mapping ([16)1 . The solution to 
the perturbation equation (JT7j) is obtained in the form (|18 )) - ()2Up assuming 
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that the parameter 9 is far from the resonance with cos# = cos 3$. Al- 
though resonant islands are important structures in symplectic maps, they 
have never been analyzed by the RG method. Here we demonstrate how our 
symplecticity-preserving RG method works near a resonant island. A similar 
procedure can be performed near other ones. 
Let us expand 9 near the resonance 

7T 



9 



+ e 0(D +e 20(2) + .. 



Eq. (|TB|) can then be expressed in the form 



20« 



x 



2e 2 9^x n , 



L n/2 x n = e[2J{x n f 
where L n / 2 x n is defined by 

L v/2 x n = x n+1 +x n ~\ 
The perturbation expansion yields 

^/2X ( ° )n = 0, 

L n/2 x (1)n = 2 J(x (0)n ) 3 - 29^x {0)n , 

L n/2 x^ n = QJ{x^ n fx^ n - 29^x^ n - 29^x^ n . 

These Eqns can be solved, yielding the result 



x 



x 



(0) n 

(1) n 



Ai n + c.c. 



-i)i n n 



J{A* 3 + 3\A\ 2 A) 



c.c. 



X 



(2)n 



zV 



J' z [ -2\A\*A+\A\ 2 A* 3 + A 



+J9 {l) [3\A\ 2 A- A* 3 ) - 



-A 



+ i n ni9 {2) A 



c.c. 



Here A(e C) is the integration constant. 

As before, we define the renormalization transformation by 



.4'' 



A + e(-i)n 
3 



J(A* 3 + 3\A\ 2 A) -# (1) aJ 
+e V|^J 2 (-2|A| 4 A + \A\ 2 A* 3 + A 5 ) + J9 {l) (3\A\ 2 A - A 



A| + £ 2 m6» (2) A + c.c. 
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Taking into account the expression 

A = A n + ein{j({A* n f + 3\A n \ 2 A n ^j - # (1) A n j, 

which relates the amplitude A to the renormalization variable A n , we obtain 
the naive RG map 

A r l +1 = Al + e(AJ(Alf -d {l) A™ 

+e 2 | - 2AJ 2 (A n l f(A n 2 f + 2 J0« [{Alf + 3A"(A™) 2 

fl(l) 2 ~i 

A n_Q{2) A n\ ^ 



2 

= A% + e[-4J(A%) 3 + 9 {1) A^ 

+e 2 | - 24J 2 (A") 2 (A™) 3 + 2J6^ ( (A™) 3 + 3(A") 2 A' 2 



-Al-e^A^y (29) 

Here A n = A™ + , and A", are real variables. can be determined 
by taking the limit e — > 



2 



dt 2 <9A 



H^{Ai,A 2 ) = (JA\-e { - l) Al/2) + (JAt-e { ' 1) Al/2). 

The smplecticity-preserving RG procedure gives the expression for H™. 
Identifying A n+1 in pH jl -(^ jl as 

A(t) + eC H A(t) + £ ^f-A{t) 

= A(t) + e{A(t), H®} + e 2 ({A(t), H^} + ±{{A{t), H^}, , 
then the trancated Liouville operator relations are 
A x (t + e) = A 1 {t)+e(AJA i 2 -6^A 2 
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2AJ 2 A\A\ + 2 J9 {1) \A\ + 3AiA 



9(1)2 a n(2) A ^ (2) 

At - 9 {2 >A 2 + 



A 2 (t + e) 



2 " -i " ' &4 2 
A 2 + e[-UAl + eWAi] 



(30) 



- 2AJ 2 A\A\ + 2 J9 {1) (A\ + 3A 2 A 



^ 1)2 „ am A ^ (2) i 

A 2 - 9 W A 1 — — 

2 2 dA x 



(31) 



Compring (|3T)|) - (f3~TJ) with ()28)) - (f2Tjj) . we have the Hamiltonian whose trajec- 
tory can approximatly interpolate the trajectory of the naive RG map, 



H = + eH (2) = J(A\ + A 



o (Al + A 2 2 ). 



This Hamiltonian system is integrable because of 1 degree of freedom. How- 
ever the solution may be complicated to write the analitic expression, which 
yields the difficulty of finding a symplecticity-preserving RG map. Instead 
of using the analytical expression to the Hamiltonian flow, we use a sym- 
plectic integrator, which is known to be a discretization method designed for 
preserving the symplecticity for Hamiltonian flows |13j. 

We split the Hamiltonian as H\ + H 2 , not + eH^ 2 \ so that we use a 
symplectic integrator, 



H = H 1 + H 2 



fl'(i) 

JA\ - —A\ ) + (JA 



y 0'(i) = 0(i) +£ 0(2). 



2 V V 2 
The Hamiltoinian Hi provides the mapping 

e rD "i : A 1 {t + r) = A 1 (t), A 2 (t + r) = A 2 {t) + (- 4J A\{t) + 0'« Ai(*))t. 
Here r is a real number. Similarly, if 2 provides 

A x (* + r) = + (^4J^(t) - fl'W^^T, A 2 (t + r) = A 2 (t). 



Therefore we can obtain the reduced symplectic map 



A'l 



n+l 



A\ + AJe 



[ A 2 + |(- 4J ^f + 0' (1) 4l)}' 
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AT 1 



-#' (1) jv4™ + E -{-UAf + #' (1) a™)} 
Al + e(^-AJAf + e' {l) A^ 

_ 4JA n+l 3 + d '(l) A n+l\ 



(32) 



(33) 



Here we have taken the following symplectic integrator: 
exp(e 



Dhi " exp(eD H2 ) exp(e^) + 0(e 3 ) 



3.3 



2 ' Ly 112/ " v 2 
Numerical results in case of a resonance 



In this subsection we present illustrative numerical results for our symplecticity- 
preserving RG method in case of near the resonance with cos 9 = cos 39. To 
show that our RG method can successfully be applied to study a resonant 
structure for the 2-dimensional symplectic map (JTHJ), we use both the reduced 
maps @2h© and flZZD up to 0(e). 




-4-3-2-101234 





-3-2-101234 



-4-3-2-101234 



Figure 1: Phase portraits of the 2-dimensional symplectic map model, with 
the parameters are e = 0.01, J = 1.0, and 9^ = 1.0: (a) the original map 
[Eq. (fTB|)]. (b) the Liouville operator approach to the RG method [Eq. (pj2~j) - 
(1331 up to 0(e) ], (c) the exponentiated RG method [Eq.flZZD up to 0(e)). 



In figures 1-2 the phase portraits to the map near the resonance with 
cos 9 = cos 36*, are depicted. Although the exponentiated RG map agrees 
well with the exact numerical result for near the origin in the phase space, it 
derivates considerably from the exact result near the resonance point. In con- 
trast, the map obtained by the Liouville operater approach to a symplecticity- 
preserving RG method can globally give a good approximation even in the 
case of the 1st RG approximation. Detailed studies of such resonant island 
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Figure 2: Phase portraits of the 2-dimensional symplectic map model, with 
the parameters are s = 0.01, J = 1.0, and 6^ = 10.0: (a) the original map 
[Eq. (US))], (b) the Liouville operator approach to the RG method [Eq. 
up to 0{e) ], (c) the exponentiated RG method [Eo.([27jl up to 0{e)\. 



structure by means of the present RG method will be published elsewhere 
for other maps including the Henon map. 

4 Conclusions 

We present the Liouville operator approach to the RG method to preserve 
symplectic structures in RG maps near elliptic fixed points of symplectic dis- 
crete systems. The symplecticity-preserving procedure is accomplished by 
comparing naive RG maps with the Liouville operator relation order by or- 
der in the small parameter and gives symplectic maps, which successfully de- 
scribes the long-time asymptotic behavior of the original systems. Although 
the exponentiation procedure has not been able to be applied to study res- 
onant islands, the Liouville operator approach to symplecticity preservation 
has given correctly reduced maps. Furthermore the advantages of this new 
method are not only a time-step parameter is not needed, but also reduced 
maps may be able to be explicit schemes. 

It is easy to see that the present symplecticity-preserving method is also 
applicable to general weakly-nonlinear symplectic maps. Other symplecticity- 
preserving RG maps are to be studied in future. 
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